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Let £' denote the set of non-normalized two-qutrit entangled states of rank four having positive 
partial transpose (PPT). We show that the set of SLOCC equivalence classes of states in £' , equipped 
with the quotient topology, is homeomorphic to the quotient R/A$ of the open rectangular box 
, R C R 4 by an action of the alternating group A5. We construct an explicit map uj : O — > £', 

where is the open positive orthant in R 4 , whose image lu(Q) meets every SLOCC equivalence 
f"^ ' class E C £'. Although the intersection uj(Q) n E is not necessarily a singleton set, it is always a 

£Sj , finite set of cardinality at most 60. By abuse of language, we say that any state in w(fi) n E is a 

canonical form of any p 6 E. In particular, we show that all checkerboard PPT entangled states 
can be parametrized up to SLOCC equivalence by only two real parameters. We also summarize 
the known results on two-qutrit extreme PPT states and edge states, and examine which other 
interesting properties they may have. Thus we find the first examples of extreme PPT states whose 
ON ' rank is different from the rank of its partial transpose. 

<N . 
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INTRODUCTION 



£ 

Q" 1 Entanglement plays the essential role in quantum information processing, such as quantum teleportation [2j, com- 
puting [32j and cryptography [l7| . Entanglement also reveals a fundamental difference between the quantum and 
^ classical world which may be detected by Bell inequallity and other means [||. In spite of these striking facts, it is a 

\ hard problem to decide whether a given quantum state is entangled jl8j . 
£SJ ■ A non-entangled state, also known as a separable state, is by definition a convex sum of product states (37j . For a 
bipartite state p acting on the Hilbert space H := 'Ha®'tIb; the partial transpose computed in an orthonormal (o.n.) 
basis {|aj)} of system A, is defined by p T = |aj)(<Jj| ® (a,i\p\aj). The dimensions of Ha and T~Lb are denoted by 
M and N, respectively. (We assume that M, N > 2.) We say that p is a k x I state if its local ranks are k and I, i.e., 
rank pa = k and rankps = I. We say that p is a PPT [NPT] state if p r > [p T has at least one negative eigenvalue]. 
Evidently, a separable state must be PPT. The converse is true only if MN < 6 [24|, [33|. The first examples of 
two-qutrit PPT entangled states (PPTES) were constructed (in purely mathematical context) by Choi and St0rmer 
in the 1980s [H, HB|. They were introduced in 1997 [25| into quantum information theory. The Choi's example has 
been generalized in [l9j . The PPTES can be used for many tasks, e.g. the extraction of distillable key [23j . 

A systematic method of constructing two-qutrit PPTES p of rank four was proposed in 1999 yf by using the 
unextendible product bases (UPB). This construction is indeed universal: it was proved in 2011 [9|, [35( that any 
such p can be constructed by using a UPB. Any such p can be converted by stochastic local operations and classical 
communications (SLOCC) into a canonical form depending on for positive parameters, see Eq. ([7]) and the subsequent 
paragraph. This canonical form is a minor modification of the one constructed in [3, Eq. (108)]. We demonstrate 



our result by using the well-known Pyramid and Tiles UPB [14[. The set of SLOCC equivalence classes of two-qutrit 
PPTES carries the natural quotient topology. We show that this topological space is homeomorphic to the quotient 
of the open rectangular box R C R 4 by an action of the alternating group A5. The checkerboard family of PPTES is 
an early well-known example of an infinite family of two-qutrit PPTES of rank four Q. We show that, under SLOCC, 
the checkerboard PPTES can be parametrized by only two real parameters. These results clarify the structure of 
3x3 PPTES of rank four. 

We say that a PPT state p is extreme if it is an extreme point of the compact convex set of all normalized PPT states. 
Every PPT state is a convex linear combination of extreme states. For convenience we say that a non-normalized 
state is extreme if its normalization is extreme. 
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We say that a PPT state p is an edge state if there is no product vector \a,b) € Ti{p) such that |a*,6) € TZ{p T ). 
Thus an edge state must be a PPTES. It is easy to see that any entangled extreme state is an edge state, but the 
converse is false. It has been proved that any 3x3 PPTES of rank four is an extreme state @. They are PPTES 
of lowest rank. Both extreme and edge states in higher dimensions have been extensively studied in recent years 
@, EL US HH • Among them, 3x3 extreme and edge states have the simplest structure. We will give a summary of 
the recent progress on these states, and point out some unknown cases. We explicitly construct states for these cases, 
including extreme PPTES of birank (5,6) and (5,7), as well as non-extreme edge states of birank (6,6), (5,6) and 
(5, 7). (We refer to the ordered pair (rank p, rank p r ) as the birank of p.) 

The content of our paper is as follows. In Sec. [XT] we provide the technical tools used in the paper. In Sec. IIIII 
we consider the set £' of 3 x 3 PPTES of rank four, which is invariant under the natural action of the group G of 
local invertible transformations. In Theorem [7] we prove that any p 6 £' can be converted under SLOCC into the 
canonical form. We also prove, see Theorem[5J that the quotient spaces £' /G and R/A5 mentioned in the abstract are 
homeomorphic. In Sec. IIV1 we construct a family of checkerboard PPTES depending on two positive parameters u 
and v, see Proposition [TU We show that any checkerboard PPTES is SLOCC equivalent to a member of this family. 
By using invariants, we find a simple characterization of checkerboard PPTES (see Proposition [T5|) . In Sec. fVl we 
summarize the recent progress on 3 x 3 extreme and edge states. We examine some of the known states of this type 
to find out what other properties they may have. E.g. we observe that there exist extreme states of birank (r, s) with 
r s. Apparently, this observation is new. The results are presented in Table I at the end of this section. We also 
propose several open problems. 



II. PRELIMINARIES 



We shall write Ik for the identity k x k matrix. We denote by IZ(p) and kerp the range and kernel of a linear 
map p, respectively From now on, unless stated otherwise, the states will not be normalized. We shall denote by 
{\i)a '■ i = 0, . . • , M — 1} and {\j)b '■ j = 0, . . . , N — 1} o.n. bases of T-La and Hb, respectively. The subscripts A and 
B will be often omitted. Any state p of rank r can be represented as 

M-l 

p=52\m®ctc,, a) 

h3=0 

where the C% are R x TV matrices and R is an arbitrary integer > r. In particular, one can take R = r. Wc shall often 
consider p as a block matrix p = C^C = [CjCj], where C = [Co C\ ■ ■ ■ Cm-i] is an R x MN matrix. Thus C\Cj is 
the matrix of the linear operator (i\Ap\j)A acting on %b- 

In physics, the density matrix p describes the systems AB. In many cases we only need to describe one system by 
the reduced density matrices pa = Trg p and ps = Tr^ p. For these matrices, using Eq. (JXJ) we have the formulae 

M-l 

Pb = E C l C ^ P A = t Tr C i C ^ i,j = 0,...,M — 1. (2) 

i=0 

It is easy to verify that the range of p is the column space of the matrix and that 

{M-l M-l \ 

E !*'> ® lw> : E °i\vi) = [ ■ ( 3 ) 
i=0 1=0 J 

In particular, if C;|j) = for some i and j then \i,j) £ kerp. 
For any bipartite state p we have 

(p r ) B = Tr A (p r ) =Tv AP = PB, (4) 
(p r ) A = Tr B (p r ) = (Tr B pf = (p A f. (5) 
(The exponent T denotes transposition.) Consequently, 

rank (p T ) A B = rankp^.s- (6) 

If p is an M x TV PPT state, then p r is too. If p is a PPTES so is p r , but they may have different ranks. 

Let us now recall some basic results from quantum information regarding the separability and PPT properties of 
bipartite states. Let us start with the basic definition. 
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Definition 1 We say that two bipartite states p and a are equivalent under stochastic local operations and classical 
communications ('SLOCC-equivalcnt, or just equivalent ) if there exists an invertible local operator (ILO) V = Va®Vb 
such that VpV^ = a 

This is the physical realization of ILO. Thus the equivalence classes of states are just the orbits under the action of 
the group G :— GL3(C) x GLa(C). It is easy to see that any ILO transforms PPT, entangled, or separable state into 
the same kind of state. We shall often use ILOs to simplify the density matrices of states. 

Let us recall from [H, Theorem 22] and [H, Theorems 17,22] the main facts about the 3x3 PPT states of rank four. 
Let M = N = 3 and let U denote the set of unextendible product bases in T-L. For {ip} £ U we denote by LT{?/;} the 
normalized state (1/4)P, where P is the orthogonal projector onto {ip} ± - We say that a subspace of H is completely 
entangled (CES) if it contains no product vectors. (We require product vectors to be nonzero.) For counting purposes 
we do not distinguish product vectors which are scalar multiples of each other. 

We give a formal definition of the term "general position" [111, Definition 7] . 

Definition 2 We say that a family of product vectors {\ipi) = \4>i) ® \Xi) ■ i £ 1} is in general position (in H) if for 
any J C I with \ J\ < M the vectors \4>j) , j £ J, are linearly independent and for any K C I with \K\ < N the vectors 
\Xk) } k £ K, are linearly independent. 

Theorem 3 (M = N = 3) For a 3 x 3 PPT state p of rank four, the following assertions hold. 

(i) p is entangled if and only iflZ(p) is a CES. 

(ii) If p is separable, then it is either the sum of four pure product states or the sum of a pure product state and a 
2x2 separable state of rank three. 

( in ) If p is entangled, then 

(a) p is extreme; 

(b) rankp r = 4; 

(c) p = V LT{-0} V* for some V £ G and some {ip} £ U; 

(d) kerp contains exactly 6 product vectors, and these vectors are in general position. 
The first assertion of the following proposition is 0, Theorem 23]. 

Proposition 4 (7\/ = N = 3) Any p £ £' is SLOCC equivalent to one which is invariant under partial transpose, 
i.e., there exist A,E> £ GL3(C) such that a := A® B p A^ (g) satisfies the equality a r = a. Moreover, we may 
assume that a = C^C where C = [Co C\ C2] and 
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The weaker form of the second assertion, with a being just real and nonzero, follows from the proof of [9|, Theorem 
23] and [l(| ■ The stronger claim that (like b, c, d) a can also be chosen to be positive will be proved in Theorem [7] 
For convenience, we say that any state a = C^C, where C is defined as in this proposition with a, b,c,d> 0, is in the 
canonical form. 

As an application, let us recall the following result 0, Theorem 24]. 

Proposition 5 (AI = N = 3) // the normalized states p and p' are 3x3 PPTES of rank four with the same range, 
then p = p' . 



III. EQUIVALENCE CLASSES OF 3 x 3 PPTES OF RANK FOUR 

The main objective of this section is to show that the set of equivalence classes of 3 x 3 PPTES of rank four, 
equipped with the quotient topology, is homeomorphic to the quotient R/A§ where R C R is a product of four open 
intervals and A§ is the alternating group of order 60 acting on R by rational transformations. We also formulate a 
test for checking the equivalence of two 3x3 PPTES of rank four. 

Recall that £' is the set of 3 x 3 PPTES of rank four, and that £' is G-invariant. The quotient space £'/G 
parametrizes the set of equivalence classes of 3 x 3 PPTES of rank four. We equip £' /G with the quotient topology 
and let 7r : £' — > £' /G be the projection map. 
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Let fi = {(a,b,c,d) : a,b,c,d > 0} be the positive orthant of R 4 . Define the map uj : — > £' by the formula 
a; (a, b, c, d) = C^C, where C = [C C\ C 2 ] and the blocks Cj are 4x3 matrices given byEq. (J7|). Our first objective is 
to prove that the map ttlo : SI — > £'/G is onto. This was proved in our previous paper Q apart from the fact that the 
sign of the parameter a was left ambiguous. The proof uses the G-invariants of quintuples of product vectors which 
we now briefly recall. 

Let (|o:i)) 4 = Q be an (ordered) quintuple of vectors in T~La such that any three of them are linearly independent. To 
any such quintuple we assign three invariants (Jf 4 , j£, J^)- These are certain complex numbers, different from and 
1, subject to the relation J^J^J^ = 1. They arc defined by the formulae 

T A _ A2,Q,4Aq,i,3 a _ A 0i i ; 4Ai j2 ,3 T A _ Al,2,4A 2;0)3 , . 

1 A A ' 2 A A ' 3 A A ' \°) 

^2,0,3^0,1,4 ^0,1,3^1,2,4 ^1,2,3^2,0,4 

where A^j.^ = det[|«i) \aj) \ctk)]- If two quintuples, say (|c*i)) and (|<a^)), have the same invariants then there exists 
an invertible linear operator Va on Ha such that VA|«i) cx |a£) for each index i. (The proportionality constants may 
depend on the index.) The converse is also valid. 

To any quintuple of product vectors (\4>i) = \a.i, /3i))f =0 in H, which are in general position, we assign six invariants 
( Ji, J2, J3, Ji , J® ■> where the (J± ,Js) are the invariants of the quintuple (|«i)) and ( J-f, J^, J-f) are those 
of the quintuple (|A))- If t w0 quintuples of product vectors, say (|0i)) and (\tpi)), have the same invariants then there 
exists an ILO, say V, such that V\(f>i) tx \ipi) for each index i, and the converse holds. 

For p € £' , ker p contains exactly six product vectors, say \4>i), i — 0, . . . , 5. Moreover, these six product vectors 
are in general position, see 0, Theorem 22]. There are in total 720 different sextuples (\<t> v i)), where n £ Sq is a 
permutation of the set {0, 1, . . . , 5}. We define the six invariants of such sextuple to be the invariants of the quintuple 
which is obtained from the sextuple by truncation, i.e., by omitting the last product vector l^s). It follows from 
0, Theorem 25] that, for all sextuples (\4>Tri)), all six invariants are real. Hence, each of the invariants belongs to 
one of the open intervals p = (0,1), P = (f,+oo) and N = (— oo,0). To each of the sextuples we associate a six 
letter symbol, with letters chosen from the set {p, P, N}. The symbol is constructed from the sequence of invariants 
( Jj 4 , J2 , J3, J\ , J2 > J3) by replacing each number by the letter designating the open interval p, P, N containing it. 
In the generic case, the 720 sextuples of product vectors will have pairwise different sextuples of invariants. However, 
only 12 different symbols arise. We refer to these particular symbols as the UPB symbols (see [|1 Table I]. Each of 
the 12 UPB symbols arises exactly 60 times. For convenience we single out one of the 12 UPB symbols, namely the 
symbol ppPNNp. 

We assume (as we may) that the chosen sextuple (|<^>i)) is of type ppPNNp. The set of all permutations tt G Sq for 
which (\4>Tri)) is of type ppPNNp is a subgroup of 5*6 isomorphic to the alternating group A5 of order 60. We shall 
refer to this subgroup as the stabilizer of the symbol ppPNNp. One defines similarly the stabilizers for other UPB 
symbols. All these stabilizers are isomorphic to A§ and are pairwise distinct. For the symbol ppPNNp, the stabilizer 
is generated by the 5-cycle a = (0, 1, 2, 3, 5) and the involution /? = (1, 2)(4, 5). Note that it permutes transitively the 
six product vectors \(j>i). From now on we set A 5 = (a, 0). 

Since J^J^J^ = 1 and Jf ' ' J§ = 1, it suffices to work with the following four invariants Jf 4 ,3% ~,J£ \Ji '• Fb r 
each 7r £ A 5 , the sextuple (|</>th)) has symbol ppPNNp and so the quadruple of its invariants ( Jj 4 , J2, j£* , J3) belongs 
to the open infinite rectangular box R = pxpxNxp(& product of four open intervals in the Euclidean space 
R 4 ). We shall see that R parametrizes the set of equivalence classes of 3 x 3 PPTES p of rank four. However, this 
paramctrization is not one-to-one; generically it is sixty-to-one. The action of A$ on sextuples of product vectors 
in ker p induces an action of A$ on R, and the equivalence classes are represented by the orbits of A§ in R. Let us 
describe explicitly the action of a and /? on the quadruple of invariants (a, 6, c, d) := ( Jj 4 , J-f, J^, J.f): 

(nhrA^ ( (i-dXfr-c) u (l-c)(l-ttcrf) 1 (l-b)(l-acd) (1 - a)(l - abd) \ 
• {a,o,c,a)^ ^ 1 _ c ^ 1 _ aM y °- (i_ cd )(fc_ c )' d (i_o)(6-c) ' (1 - 06) (1 - acd) J ' {1 

. / 1 - c (1 - ab) - ce b (1 - ab) - ce (1 - a)(l - abd) \ 

P . (aAM)^-—, (1 _ c)(1 _ flM) , c'(l-6)(l-aW)' °' (I - ab) - ce J' (10) 

where e = (1 — a) + ad(l — b). The meaning of, say, the first formula is that its right hand member is the quadruple 
of invariants of the sextuple (\(f> ai )) = (\(j> 5 ), |0 O ), |0i), ^2), \<f>4,), \<h))- 

From the preceding discussion we obtain the following test for equivalence. 

Theorem 6 Let p,p' G £' and let (\4>i)) and (\4>'j)) be sextuples of product vectors in kerp and ker//, respectively. 
Then p and p' are SLOCC equivalent if and only if there exists a permutation 7r € Sq such that the sextuples {\<j>m)) 
and {\4>'i)) have the same invariants. 
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Proof. Necessity. Let p' = V pV\ where V is an ILO. Since p\<fti) = and kcrp = W kcrp', we have V^l^) = Ci\(j>„i) 
for some permutation 7r and some nonzero scalars Cj. It follows that the sextuples (|<^ W j)) and (|</>i)) have the same 
invariants. 

Sufficiency. Since (\4>ni)) and (\4>'i)) have the same invariants, there exists an ILO, say V, such that V"f|(^) = Ci\4>iri) 
for some permutation ir and some nonzero scalars Cj. Since the span ker// and the |</>j) span ker/?, we have 
ker p' = kcrp = kei(VpV^). It follows that the states p and VpV^ have the same kernel and range. By 
Proposition [5] these two states are equivalent, and so are p and p'. □ 

We shall now prove that each point of R is the quadruple of invariants of some p G £' and that we can indeed 
require in Proposition [4] that a > 0. In view of Theorem [6] this implies that the map ttlo : O — > £' /G is onto. 

Theorem 7 Every a G £' is SLOCC equivalent to a state p = u(a,b,c,d) (where all a, 6, c, d > 0J. For eac/i point 
(x, y, z,w) £ R there exists a state /)££' suc/i t/ia£ (x, y, z, u>) = ( J-j 4 , J^ 1 , Jf, J.f) for some sextuple of product vectors 
in ker p. 

Proof. To prove the first assertion we choose a sextuple of product vectors (\4>i))i=o m ker a with invariants of type 
ppPNNp. Thus we have (Jj 4 , J A , J 2 7J3) £ -R- Let us define the positive numbers b, c, d by 



(l-J^Xi-Ji 4 ^) 



j a (i-jfj£)(i-jfj?jBy 
1 {i-J?J§){j£-J§Ji) 



c 



Jf (1-J 3 B )(1-W) 



(11) 

(12) 



, (1-J 2 fl)(l-J 3 *) 



(It is easy to see that the right hand sides are positive.) 
We can now define a > by the formula 



bed (l - J*)(i - J A J 2 B J 3 B )(J 2 A - JfJ£) 



Ji (1- J A )(1- J^jf)2 



(14) 



We claim that p := w(a, 6, c, d) and cr are equivalent. It suffices to show that we can choose a sextuple of product 
vectors in ker/? having the same invariants as the sextuple (\4>i))- We observe that kerp contains the product vectors 
ii) for i = 0, 1, 2. We have to find the remaining three product vectors in this kernel. We proceed as in the proof of 
Theorem 25]. We introduce the cubic polynomial 

f(z) = abz(cz-l-d 2 ){c-(l + c 2 )z) + d(cz-l)(b 2 c-(l + b 2 + b 2 c 2 )z). (15) 

It has three real roots Z\,Z2, 23 such that 

z 3 <0, A < zx < c/(l + c 2 ), l/c< z 2 < (l + d 2 )/c, (16) 

where A = 6 2 c/(l + b 2 + b 2 c 2 ). Explicitly, they are given by the formulae 



J_i_ 1 J A Jj 

c '(i-jfjfjfy 



2 i = — "7i tA tB T By » ( 17 ) 



_ 1 1 - Jj 4 Jf 

22 ^ c ' (l-jfjfjfy ( 8) 

I . (i-^)d-j^f) (19) 

[The verification that these are indeed the roots of f(z) is a very tedious job if done by hand, but it is trivial for a 
package for symbolic algebraic computations such as Maple.] 

The three additional product vectors in the kernel of p are \ipi)(i = 1,2,3) defined by 0, Eq. (110)]. By using 
the above expressions for a,b,c,d and 21,22,2:3 and the formulae [9j, Eqs. (119-120)], we can compute the invariants 
of the sextuple (|00), |11), |22), IV^), \i/>3))- We obtain that they are the same as the invariants of the sextuple 
(\4>i))- This proves our claim and completes the proof of the first assertion. 

To prove the second assertion, let us define a smooth map $ : R — >• f2. We set <&(x, y, z, w) = (a,b,c,d), where 
a, b, c, d are defined by the formulae (jTTJ) - (jT^J) in which we replace J-j 4 , J A , J 2 , j£ with x, y, z, w, respectively. The proof 
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of the first assertion shows that the quadruple of invariants of the state p, constructed above, is exactly (x, y, z, w). 
Hence, the second assertion is proved. □ 
It is not hard to show that the above map $ is a local diffeomorphism. We believe that it is a (global) diffcomorphism, 
but this question is beyond the scope of this paper. However, let us observe that if (x, y, z, w) € $ _1 (a, b, c, d), then 
(x, y, z, w) is a quadruple of invariants of some ordered sextuple of the six product vectors in the kernel of w(a, b, c, d). 
Moreover, the type of this sextuple is ppPNNp. Consequently, the fibre $ -1 (a, 6, c, d) is finite of cardinality at most 
60. 

Corollary 8 The set, £, of all normalized 3x3 PPTES of rank four on T-L has dimension 36. 

Proof. This set is a real semialgebraic set and so it has a well defined dimension, see e.g. @. For convenience 
we shall work with non-normalized states, and so the dimension will increase by one. More precisely, we have 
£' = {Xp : A > 0,p G £}. Denote by X the image of u>, X = {uj(a,b,c,d) : (a,b,c,d) £ fl} C £'. Each equivalence 
class E C £' is an orbit of G. It follows from [j| Proposition 27] that the identity component of the stabilizer of p € X 
is the 3-dimensional subgroup {(2:1/3,22/3) : \ziz 2 \ = 1} of G. Hence, DimE = DimG — 3 = 33. Each equivalence 
class E intersects X in at least one and at most 60 points. Consequently, Dim£' = Dim A + DimE = 37. The 
normalization decreases the dimension by 1, thus Dim£ = 36. □ 
We point out that this agrees with the computation performed in [22j . 

Let p : R — > R/A$ be the canonical projection map onto the quotient of R by the action of A5. As usual, the quotient 
space R/A5 is equipped with the quotient topology. We cannot define a map £' — > R because there is no natural 
choice of ordering the six product vectors in the kernel of a state p £ £' . However, this ambiguity disappears when we 
replace R by R/A$ and specify that the sextuple of product vectors in kerp is chosen to be of type ppPNNp. Thus 
we obtain a continuous map q : £' — > R/A5. We now use the universal property of the quotient maps tt : £' — > £'/G 
and p : R — > R/A$ to prove the following. 

Theorem 9 The quotient spaces £' /G and R/A§ are homeomorphic. 

Proof. Since the map q : £' — > R/A$ is continuous and is constant on the fibers of the map tt : £' — > £'/G, there 
exists a unique continuous map / : £' /G — > R/A§ such that q = fix. At the end of the proof of the previous theorem, 
we have introduced the smooth map $ : R — > 0. Since the map 7rw$ : R — > £' /G is continuous and constant on the 
fibers of the map p : R — > R/A$, there exists a unique continuos map g : R/A5 — > £' /G such that 7tcj$ = gp. The 
universal property also implies that fg and gf are the identity maps. Hence, / is a homeomorphism. □ 
We shall now examine the two classical examples of UPB, namely Pyramid and Tiles. 

Example 10 The quadruple of invariants for the Pyramid UPB is (-(y/E-l)/2, (3 + y/E)/2, (3-a/5)/2, (y/5 + l)/2) 
with symbol NPNPpP. If we rearrange the five vectors of this UPB in order [2, 1, 4, 3, 5], the quadruple of invariants 
becomes ({V~E- l)/2, (V5-l)/2,-(VE+l)/2, (3- \/5)/2). So it has the desired symbol ppPNNp and the point with 
these coordinates is in the region R. 

We claim that this point is the unique fixed point of the transformation a as given by Eq. Indeed, by equating 
the two members of Eq. (|9|), we obtain a system of four equations. The first two equations can be solved easily and 
give 

l-b-ab + a 2 b 2 a + b-1 

° = (l-a)(l-a-by l-b{l-aby { ' 

Then the remaining two equations factorize, and after droping the factors that cannot vanish when a, b € (0, 1), we 
obtain the system of two equations: 

a 2 b 3 - ab 2 + ab - a - b+ 1 = 0, (21) 
a 2 b 3 - a 2 b 2 - a 2 b - 2ab 2 + Sab + b - 1 = 0. (22) 

The resultant of these two polynomials with respect to the variable a is 6(1 + b 2 )(\ — b) 3 (l — b — b 2 ). It has a unique 
root in the interval (0, 1) namely b = (v5 — l)/2. The proof of our claim can now be completed easily. 

Moreover, one can verify that the transformation fj fixes the unique fixed point of a. Hence, this point is fixed by 
all elements of A5. Note that this observation agrees with the fact mentioned in Q that the stabilizer of Pyramid 
PPTES in the group PGL 3 x PGL 3 is the alternating group A 5 . □ 



Example 11 Let us now consider the Tiles UPB. The quadruple of its invariants is (—1/2,2,1/3,3/2) again 
with symbol NPNPpP. By rearranging the five vectors in order [2,1,4,3,5], we obtain the quadruple of invariants 
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(1/2, 2/3, -1, 1/2) having the desired symbol ppPNNp. The orbit of A 5 containing the point (1/2, 2/3, -1, 1/2) G R 
has size five. The other four points of this orbit are (1/2,1/2,-2,1/4), (2/3,1/2,-2,1/2), (2/3,3/4,-3,1/3) and 
(3/4, 2/3, —1, 1/3). Each of these five points (x, y, z, w) gives the parameters (a, b, c, d) = $(x, y, z,w) € Q such that 
the state oj(a, b, c, d) € £' is equivalent to the PPTES obtained from Tiles. For instance, by substituting the coordi- 
nates of the first point for the invariants J-f 4 , J® , J® in Eqs. fjl 141 14[) . we obtain that a = 7^21/27, b = 2/3a/5, 
c = v / 2T/2 and d = 2\/5. □ 

At the referee's suggestion, we shall also consider the one-parameter generalization of Choi's PPTES. 
Example 12 This one-parameter family p\ of 3 x 3 PPTES of rank four is given in [lj| p. 169] by its density matrix 



1 



A 2 
1 



A 2 
1 



1 



1 



1 



(23) 



(The blank entries are zeros.) The parameter A is assumed to be positive and different from 1. For simplicity, we 
shall assume that < A < 1. The six product vectors in kcrp,\ are the tensor products of the corresponding columns 
in the following two 3x6 matrices 



1 





A 


1 


A 


1 





-A 





A 


1 








1 

-A 



-A 
1 





-A 
1 



1 A 1 
1 A 
-A 1 A 



(24) 



This sextuple of product vectors has symbol pNNPpp. After interchanging the third and fourth columns of both 
matrices, the symbol becomes ppPNNp and a computation shows that the quadruple of invariants of the reordered 
sextuple is given by 



A 3 



1-A 3 
1 + A 3: 



1° — - 



1 + A 3 
1-A 3 



J B - 



2A 3 



1 + A 3 ' 



(25) 



By applying the formulae (fTTj) - (fl"4"]) . we conclude that p\ is equivalent to the state a defined in Proposition |H with 
a, b, c, d > given by 



2A 6 



1 + A 6 ' 



(3 + A 6 )(l + 3A 6 ) 
(1-A 6 ) 2 



d 2 = 



1 + A 6 ' 



bed (1- A 12 )(1 + 3A 6 ) 
2A 6 (3 + A 6 ) 2 ' 



(26) 

□ 



Finally we show a simple result on PPTES. 
Lemma 13 Given p 6 there is no 3 x 3 extreme state a of rank five, such that 72. (<r) 



ker p. 



Proof. Suppose such a exists. We may assume that p = uj(a,b,c,d) for some a,b,c,d > 0. Using Eq. (|7|), the 
product vector 1 00} belongs to lZ(a) and lZ(a r ). So a is not an edge state, and we have reached a contradiction. □ 



IV. CHECKERBOARD STATES 

One of the early examples of PPTES was the multi-parameter family of checkerboard states constructed in 0, 
see also [H]. We will show that, up to equivalence, all checkerboard PPTES can be parametrized by just two real 
parameters. By using invariants, we have devised a test for checking whether an arbitrary 3x3 PPTES of rank four 
is equivalent to a checkerboard state. 
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We define the checkerboard states (see [1, section 7]) to be the states p which can be written as p = C^C, where 
C = [Ci C2 C3] and the Cj are complex matrices having the following form: 




C 2 = 




Co = 




(27) 



Proposition 14 Every checkerboard PPTES is equivalent to one with 



Cx = 




Co 




c. 




u, v > 0. 



(28) 



Consequently, there exist 3x3 PPTES of rank four not equivalent to any checkerboard state. 

Proof. Let p = C^C be a 3 x 3 PPTES of rank four, where C = [Cx C 2 C 3 ] with the C t given by (J27J). Then any 
linear combination of the Ci must have rank at least two. By replacing C\ with a suitable linear combination of C\ 
and C3, we can assume that am = dj. By applying an ILO and premultiplying C by a unitary matrix, we may assume 
that 



Ci 




C 2 




c. 




(29) 



Since p v > 0, we must have e 



0. By using another ILO, we can further simplify these matrices to obtain 



C x = 




C 2 = 




C, = 




(30) 



As p is entangled, 7Z(p) must be a CES and so chlr ^ 0. By using an argument from the proof of [H, Theorem 28], 
we obtain that h = rk* , |c| = 1 and / = cr*k/k*. Thus 



Cx = 



/l 0\ 




f 


c 


10 


c 2 = 


1 











cr*k/k* 


\ / 




W) 




















rk* 




1 fc 





; ) 




V 


r 


0/ 



(31) 



Let us choose U2 such that u\ = c and premultiply C with the unitary matrix diag(l, 112, 1*3, U4), where M3 = |r|/e*/r*fc 
and U4 = U2\k\/k. Next we postmultiply each Ci with diag(l, u^/kli and multiply C2 and C3 with and 

|rfc|/rfc*, respectively. (Note that these transformations form an ILO.) Now we have 



Cx 



/ 1 

i/H 
000 
V 



Co 




c. 



/ 

\k\ 

\k\ 1 

V 1 



(32) 



Thus the first assertion is proved, and the second follows from the first because the matrices ([28]) depend only on two 
positive parameters. □ 
We can now characterize the checkerboard PPTES up to equivalence. 



Proposition 15 A 3 X 3 PPTES of rank four is equivalent to a checkerboard state if and only if, for some ordering of 
the product vectors in its kernel, its invariants are of type PNNpPP and satisfy the equations J 2 = J% and J® = J3 ■ 

Proof. Necessity. Let p be a 3 x 3 checkerboard PPTES of rank four and (J-j 4 , J^, J%, J] 8 , J^, Jjf) its invariants 
with symbol PNNpPP. By Proposition [H we have p = C^C where C = [Cx C 2 C 3 ] and the d are given by (|28p . We 
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have to find the six product vectors in ker p. First observe that the product vectors |</>o) = 1 02) and |0s) = 1 20) — v\22) 
belong to ker p. Let us exhibit the remaining four. 

Let ±Xi, ±x 2 be the roots of the biquadratic polynomial f(x) = u 2 (l + v 2 )x A — (u 2 v 2 + 2u 2 + v 2 )x 2 + u 2 . Since 
/(l) = —v 2 < 0, we may assume that < x\ < 1 < x 2 . For any real t let \(f>(t)) = \a(t)) <X> \fi(t)), where 
\a(t)) = v\0) + tv\l) + u(t 2 - 1)|2) and \(3(t)) = uvt 2 \0) - tv\\) + (t 2 - 1)|2). One can easily verify that the product 
vectors \<f>i) = \(f>(xi)), \<p 2 ) = \<p(-xi)), |03> = \<f>(x2)), \<Pa) = \<p(-x 2 )} belong to ker p. 

A computation shows that the invariants of the sextuple (|0i))f =o arc (1/ p 2 , —p, —p, 1/A 2 , A, A), where A = {x 2 + 
Xi)/(x2 — xi) > 1 and p = A(l — xix 2 )/(l + x\x 2 ) < A. As x\x\ = 1/(1 + v 2 ) < 1, we have p > 0. From the equality 
x\ = (A — /i)(A — 1)/(A + p)(X + 1) we deduce that p < 1. Hence, the invariants indeed have the required symbol 
PNNpPP and the equalities j£ = j£ and Jf = Jf hold. 

Sufficiency. Now assume that tr is a 3 x 3 PPTES of rank four and that ( Jj 4 , J^, J$, J\ \ J 2 , J-f ) are the invariants, 
having the symbol PNNpPP, of some sextuple of product vectors in ker a. We also assume that j£ = and J 2 = J3 , 
and so we can write these invariants as (1/p 2 , —fi, —p, 1/A 2 , A, A), where A > 1 and p > 0. Since the first letter of the 
above symbol is P, we have p < 1. We set 



It is not hard to verify that 



Xu X — a A" + 1 . 
>0, w = 2-- ■ . (33) 



X — p ' X + p A 2 — 1 



' 1+ " 2 »'-- 1 '- 1 ° ^ ! w^ > - < 34 > 

and so we can define u > by the equation u 2 ((l + v 2 )(w — 1) — 1) = v 2 . 

Let p = C^C be the checkerboard state where C = [C\ G 2 C3] and the Cj are given by (|28[) with the parameters 
u and w as defined above. The biquadratic polynomial f(x) = u 2 (l + v 2 )x 4 — (u 2 v 2 + 2u 2 + v 2 )x 2 + u 2 has roots 
±xi,±a:2) where x\,x 2 > are given by 



(A-/j)(A-l) 2= (X-p)(X + l) 

(x+p)(x + iy X2 (x+p)(x-iy 



(35) 



The computation performed in the first part of the proof shows that (1/p 2 , —p, —p, 1/A 2 , A, A) are also the invariants 
of a sextuple of product vectors in ker p. Hence, p and a are equivalent. □ 
We remark that if we use the invariants of type ppPNNp (instead of PNNpPP) then the two equations in the 
proposition should be replaced by j£ + J% = 2J^J^ and J 2 + J3 =2J 2 . Since in Theorems 16171 we use only the 
invariants J-f 4 , J% , J 2 and J-f , we mention that the former equation can be rewritten as J 2 {1 — J\ J 2 ) = 1 — J 2 ■ 



V. EXTREME AND EDGE 3x3 STATES 



In recent years extreme and edge 3x3 states have been extensively studied. A primary reason is that PPTES do 
not exist in any bipartite space of smaller dimension (2~H. [26l |27|. On the other hand, many important problems arc 
open for 3(8)3 systems. For example, a 3 x 3 PPT state of rank four is separable if and only if there is a product vector 
in the range of this state @- But there is no analytical criterion for deciding which 3x3 states of rank larger than 
four are separable. It is also unknown whether two-qutrit entangled states of rank larger than three are distillable 
0,3- Characterizing 3x3 extreme and edge PPTES may provide a better understanding of these problems. 

We begin with extreme states. There is a simple necessary condition for extremality of states (3l| . 

Proposition 16 Let p be an M x TV PPT state of birank (r, s). If r 2 + s 2 > AI 2 N 2 + 1 then p is not extreme. 

By this proposition, the birank of an extreme 3x3 PPTES p must be (up to ordering) one of the following five 
pairs: 

(4,4), (5,5), (6,6), (5,6), (5,7). (36) 

Examples for the first three pairs have been known for some time. First, the extreme states p of birank (4, 4) were 
thoroughly analyzed and characterized by UPB 0, [35[. Second, the edge states p of biranks (5,5) and (6,6) have 
been analytically constructed [l3| . They were proved to be extreme in Ref. [2(| [23. Third, families of edge states p 
for the remaining two pairs, (5,6) and (5, 7), have been constructed very recently [30j. In Table I at the end of this 
section, we will show that for suitable parameter values, these edge states turn out to be extreme. They are the first 
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examples of extreme states p with rank p 7^ rank p T . Hence, each of the pairs (j36]) is the birank of some 3x3 extreme 
state. 

Next we consider the 3x3 edge states p of birank (r, s). 

By Q we have r = 4 if and only if s = 4. It is known that r + s < 14 [2^, [30|. (It has been claimed in [34[ that 
r + s < 13. However, it was shown in [2^, [30[ that this claim is false.) The birank (r, s) cannot be equal to (5, 9) 
because any 5-dimensional subspace of T-L contains at least one product vector. 

We conclude that, up to ordering, (r, s) must be one of the following eight pairs: 

(4,4), (5,5), (6,6), (5,6), (5,7), (5,8), (6,7), (6,8). (37) 

Recall that there exist extreme states for the first five pairs. As any extreme state of rank bigger than one is also 
an edge state, it remains to consider only the last three pairs in this list. Such edge states have been constructed 
recently in [301 ] . Since the examples for the first five pairs in (|37[) are extreme states, it is an interesting question to 
ask whether in these five cases there exist 3x3 edge states p which are not extreme. 

We have constructed such states p (see Table I) for pairs (6,6), (5,6) and (5,7). Since every PPTES of birank 
(4, 4) is extreme, only the case (5, 5) remains in doubt. We have tested a few known 3x3 edge states of birank (5,5), 
but all of them turned out to be extreme (see Table I). We conjecture that any 3x3 edge state of birank (5, 5) is 
extreme. This should be compared with the known result that any 2x4 PPTES of birank (5, 5) is extreme The 
3x3 case is more challenging because there exist 3x3 PPTES of birank (5, 5) which are not edge states. An example 
is p = a + e|00)(00| where a is a 3 x 3 PPTES of birank (4, 4), and e > is sufficiently small. 

Another interesting problem is whether our results can be extended to higher dimensions. For instance, examples of 
extreme M x N PPTES of rank M + N — 2 have been constructed recently [ll[ . One may ask how many parameters 
are necessary to describe such states. 

Table I of 3 x 3 edge and extreme states p. 
(Computations performed using Maple.) 

1. Birank (5, 5) 

(a) [H Section II] 
Extreme state. 

kerp: 2 product vectors. 

TZ(p): Infinitely many product vectors. 

7Z(p T ): Infinitely many product vectors. 

(b) [HHH with s = 1 and t = 2. 
Extreme state. 

kerp: 2 product vectors. 

IZ(p): Infinitely many product vectors. 

7Z(p T ): Infinitely many product vectors. 

(c) [il Sec. 4] with b = 2 and 6 = tt/6. 
Extreme state. 

kerp: CES. 

IZ(p): 6 product vectors. 
7£(p r ): 6 product vectors. 

2. Birank (6, 5) 

(a) [2l| with s = 2, 

Edge state, but not extreme, 
kerp: 2 product vectors. 
lZ(p T ): 3 product vectors. 

(b) [11 Sec. 4] with b = 2, 6 = tt/4 and r = 0. 
Extreme state. 

kerp: CES. 

lZ(p r ): 6 product vectors. 

3. Birank (6, 6) 

(a) O, Section III] 
Extreme state, 
kerp: CES. 
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(b) [H Sec. 4] with b = 2, 6 = tt/6, (£]??) = and (r?|C) - (CIO = 1- 
Edge state, but not extreme, 
kerp: CES. 

4. Birank (7, 5) 

(a) [H with s = 2. 

Edge state, but not extreme, 
kerp: CES. 

TZ(p T ): 6 product vectors. 

(b) HI Sec. 4] r = 3/7, 6 = 2 and 9 = arccos(9/14). 
Extreme state. 

kerp: CES. 

lZ(p T ): 6 product vectors. 

5. Birank (7, 6) 

(a) [HI with s = 2. 

Edge state, but not extreme, 
kerp: CES. 

6. Birank (8, 5) 

(a) [H Sec. 4] with b = 2, 6 = tt/4 and r = 3/7. 
Edge state, but not extreme, 
kerp: CES. 

lZ(p T ): 6 product vectors. 

7. Birank (8, 6) 

(a) 0, Sec. 3] with b = 2 and 6 = tt/4. 
Edge state, but not extreme, 
kerp: CES. 
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